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FOREWORD 


This  report  describes  the  results  of  an  experimental  program  to 
investigate  mechanisms  of  stress  corrosion  in  glasses  and  ceramics  and 
to  determine  effects  of  composition  and  microstructure  on  the  fracture 
of  piezoelectric  materials.  _  ...  _ 
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INTRODUCTION 

x 

One  of  the  main  reasons  for  strength  degradation  of  structural 
ceramic  materials  is  stress  corrosion  cracking.  Caused  mainly  by  water 
in  the  environment,  stress  corrosion  cracking  results  in  subcrltical 
crack  growth  that  leads  to  delayed  failure  of  ceramic  materials.  As  a 
consequence  of  this  crack  growth,  ceramic  components  that  are  subjected 
to  external  stresses  are  liable  to  fail  after  an  indeterminate  period  of 
time,  leading  to  a  considerable  uncertainty  in  the  safe  design  stress. 

One  of  the  objectives  of  modern  ceramics  technology  is  to  reduce  the 
uncertainty  associated  with  structural  design,  and  hence,  to  improve  our 
capabilities  of  designing  ceramic  components  that  exhibit  reliable 
mechanical  behavior.  This  objective  can  be  accomplished  in  part  through 
a  better  understanding  of  the  fundamental  processes  that  occur  during 
stress  corrosion  cracking. 

Recent  advances  in  the' area  of  structural  design  with  ceramic 
materials  have  lead  to  several  new  procedures  for  evaluating  materials 
performance  so  that  lifetime  under  load  can  be  predicted.  These  techniques 
are  all  based  on  the  science  of  fracture  mechanics,  which  provides  a 


theoretical  basis  for  understanding  how  crack  growth  leads  to  structural 


failure.  Fruitful  use  of  these  techniques  requires  that  crack  growth 
behavior  of  a  material  be  characterized  so  that  long  term  strength 
behavior  can  be  predicted.  This  characterization  can  be  accomplished 
either  by  the  use  of  fracture  mechanics  techniques  (measurement  of  crack 
velocity  as  a  function  of  stress  intensity  factor),  or  by  the  use  of 
strength  techniques  (e.g.  stress-rupture  experiments).  Current  studies 
on  the  subject  of  stress  corrosion  cracking  indicate  a  pressing  need  for 
basic  research  to  elucidate  the  chemical  and  physical  processes  that  are 
pertinent  to  the  process  of  stress  corrosion  cracking.  This  information 
could  then  be  used  to  develop  basic  models  which  would  improve  our 


i 


capability  of  making  long  term  failure  predictions  for  specific  engineering 
applications. 

Our  present  state  of  understanding  of  stress  corrosion  cracking  of 
glass  and  other  ceramics  is  primitive  in  the  sense  that  although  water 
is  known  to  he  the  main  stress  corrosion  agent,  the  actual  chemical 
reactions  causing  the  stress  corrosion  cracking  are  virtually  unknown. 
Hence,  while  ve  can  measure  crack  growth  rates,  or  strength  degradation 
rates  for  specific  ceramic  materials,  our  lack  of  knowledge  of  stress 
corrosion  mechanisms,  prevents  us  from  selecting  materials  compositions, 
or  microstructures  that  have  a  maximum  resistance  to  crack  growth.  This 
same  lack  of  knowledge  prevents  us  from  applying  fracture  mechanics 
theory  to  predict  component  lifetime,  when  the  lifetime  required  greatly 
exceeds  certain  time  limits  established  by  the  experimental  crack  growth 
or  stress  rupture  data.  The  work  reported  herein  was  obtained  as  part 
of  a  program  designed  to  further  our  understanding  of  the  stress  corrosion 


mechanism. 
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Abstract 


The  mechanical  strength  of  many  glasses  and  ceramic  materials 
decreases  with  time  under  static  loading  and  ambient  environments.  This 
strength  loss  is  associated  with  slow  growth  of  pre-existing  surface 
flaws  due  to  stress-corrosion  by  water  from  the  surrounding  environment. 
In  order  to  make  strength  predictions  for  ceramic  components  with  long 
term  applications,  it  becomes  necessary  to  understand  the  mechanism  by 
which  this  stress-corrosion  occurs.  We  have  studied  stress-corrosion  In 
vitreous  silica  exposed  to  water  and  several  non-aqueous  environments 
and  report  here  that  environments  vhich  enhance  stress-corrosion  crack 
growth  contain  active  groups  with  electron  donor  sites  on  one  end  and 
proton  donor  sites  at  the  other.  These  results  suggest  a  detailed 
chemical  model  for  the  interaction  of  the  environment  with  mechanically 
strained  bonds  in  the  solid  at  the  tip  of  a  crack.  The  proposed  model 
for  stress-corrosion  crack  growth  also  has  implications  for  the  long 
term  strength  behavior  of  a  wide  variety  of  brittle  materials. 


1  v 


Previously  published  results^-  have  shovn  that  slow  extension 
—2 

(<10  m/s)  of  a  crack  in  vitreous  silica  is  facilitated  by  water  in  the 
environment.  This  effect  can  be  demonstrated  in  a  plot  of  crack  velocity 
(V)  versus  stress  intensity  factor  (K^) ,  where  is  a  measure  of  the 
stress  field  concentrated  by  the  crack  tip.  Crack  velocities  were 
obtained  by  visual  observation  of  crack  extension  in  a  double  cantilever 
beam  type  specimen.  The  important  features  of  such  a  plot  for  vitreous 
silica  (Figure  1)  are: 

(1)  In  a  water  free  environment  such  as  vacuum,  spontaneous  failure  is 
observed  when  reaches  a  critical  value. 

(2)  When  small  amounts  of  water  vapor  are  introduced,  crack  growth  in 
the  low  velocity  regime  is  enhanced  (curve  is  shifted  to  lower  Kj 
and  V  is  a  measurable  function  of  Kj).  The  relative  amount  of 
enhancement  is  related  to  the  chemical  activity  of  water  in  the 
environment.2 

(3)  A  plateau  in  crack  velocity  is  observed  in  water  containing  environments 

(e.g.  moist  Ng).  This  plateau  has  been  linked  to  crack  growth 

limited  by  the  rate  of  water  transport  to  the  tip  of  the  moving 
2 

crack.  At  sufficiently  high  V  the  plateau  behavior  ends  and  rapid 
fracture  is  observed. 

(U)  Liquid  water  enhances  crack  growth  over  the  entire  range  of  crack 
velocity  (10-^  -  10~2). 

Several  generalized  mechanisms  have  been  proposed  to  explain  the 
effect  of  water  containing  environments  on  crack  growth  in  glass. 

Orowan  suggested  that  the  environment  could  reduce  the  energy  of  the 
fracture  generated  surface  thereby  creating  a  more  energetically  favorable 
process  for  cracking.  Charles  and  Hillig  have  discussed  stress  enhanced 
dissolution  of  material  at  the  crack  tip  and  its  effect  on  the  strength 
of  alkali  silicate  glass.  Finally,  Wiederhorn^  used  dielectric  charge 


1 


compensation  in  a  model  to  explain  environmentally  enhanced  fracture  of 
glass.  This  communication  will  describe,  on  a  molecular  scale,  a  specific 
chemical  reaction  between  strained  bonds  in  vitreous  silica  and  water 
which  can  be  used  to  explain  environmentally  enhanced  crack  growth. 
Experimental  results  in  support  of  this  model  will  also  be  presented. 

Since  we  wish  to  present  a  molecular  picture  of  the  chemical,  reaction  • 
leading  to  fracture,  a  brief  description  of  the  structure  and  bonding  in 
the  solid  and  the  environment  is  in  order. 

Network  silicates  (quartz,  coesite,  vitreous  silica,  etc.)  are 

composed  of  [SiO^]~^  tetrahedral  units  linked  at  their  corners  by  shared 

or  bridging  oxygen  atoms  (0^br ^ ) .  On  an  atomic  scale,  the  bridging  Si  - 

0  -  Si  bond  is  the  stress  apporting  member  of  the  structure  and  will  be 

important  to  the  fracture  process.  Although  some  controversy  remains 

over  the  exact  nature  of  the  bridging  bond,  recent  ab-initio  MO  calculations 

for  orthosilicic  and  pyrosilicic  acids  indicate  that  both  o-bonds  and  tt- 

3 

bonds  are  formed.  It  is  suggested  that  an  sp  hybrid  orbital  from  Si  o- 
bonds  with  an  (s  +  Apz)  hybrid  from  and  that  lone  pair  px  and  py 

orbitals  from  are  involved  in  ir-bonding  with  all  five  d  orbitals 

from  Si.  These  results  are  in  good  agreement  with  experimental  measurements 
for  ortho  and  pyrosilicic  acids  as  well  as  network  silicates. 

At  the  crack  tip  in  a  stressed  solid,  highly  concentrated  stress 
fields  are  produced.  Continuum  approximations'’  indicate  that  the  bridging 
bond  experiences  strains  >20$.  The  effect  of  this  strain  on  the  bonding 
molecular  orbitals  is  not  easy  to  predict  but  can  be  discussed  in  terms 
of  a  decrease  in  the  overlap  between  atomic  orbitals  thus  increasing 
their  availability  for  bonding  with  other  species. 


The  structure  and  bonding  in  the  water  molecule  is  established  in 
the  literature.  Oxygen  atomic  orbitals  2s,  2px  and  2py  form  three 
hybrid  orbitals  two  of  which  a-bond  with  hydrogen  atoms.  The  remaining 
lone  pair  orbitals  Cone  hybrid  and  the  remaining  ps)  are  directed  away 
from  the  hydrogen  atoms.  This  arrangement  results  in  a  positive  charge 
center  at  the  hydrogen  end  of  the  water  molecule  and  a  negative  charge 
center  on  the  opposite  end  of  the  molecule  where  the  lone  pair  orbitals 
are  located. 

Interaction  between  strained  bridging  bonds  in  silica  and  water 
from  the  environment  can  be  represented  by  a  three  step  process  (see 
Fig.  2). 

Step  1  -  A  water  molecule  from  the  environment  attaches  to  a  bridging 
Si-O-Si  bond  at  the  crack  tip.  The  water  molecule  is  aligned  by: 

1)  formation  of  a  hydrogen  bond  with  the  atom  and 

2)  interaction  of  the  lone  pair  orbitals  from  0^  with  the  Si 

atom.  The  lone  pair  orbital  interaction  may  involve  either  Van  der  Waals 

attraction  or  some  covalent  bonding  with  unoccupied  d  orbitals  of 

Si. 

Step  2  -  A  concerted  reaction  occurs  in  which  proton  transfer  to 

the  is  accomplished  simultaneously  with  electron  transfer 

from  the  0,  x  to  the  Si  atom.  As  a  result  of  this  reaction,  two 
(v) 

new  bonds  are  formed  (one  between  0^  and  Si,  and  one  between 
hydrogen  and  0^^;  the  original  briding  bond  (between  and 

Si)  is  destroyed. 

Step  3  -  Rupture  of  the  hydrogen  bond  between  0^  and  transferred 
hydrogen  occurs  to  yield  surface  Si-O-H  groups  on  each  fracture 
surface.  Since  the  hydrogen  bond  is  weak,  this  step  is  expected  to 


occur  Immediately  after  proton  transfer.  (Budd  has  proposed  a 
similar  mechanism  for  the  dissolution  of  slica  in  vater,  however, 
it  should  he  noted  that  the  present  model  for  crack  growth  is  not 
contingent  on  the  removal  of  material  from  the  fracture  surface). 

Q 

Reaction  rate  theory  can  he  used  to  show  how  the  proposed  mechanism 
for  bond  rupture  leads  to  enhanced  crack  growth  rates.  In  vacuum,  the 
entire  energy  of  the  bridging  bond  must  be  overcome  before  bond  rupture 
can  occur.  (This  energy  is  supplied  in  part  by  mechanical  strain  as 
well  as  by  thermal  contributions).  When  water  is  present,  the  height  of 
the  energy  barrier  to  bond  rupture  is  reduced  because  charge  transfer 
(Step  2)  provides  an  alternate,  lower  energy,  reaction  path.  Accordingly, 
the  activated  complex  is  that  structure  which  exists  during  the  charge 
transfer  reaction.  Since  the  rate  of  a  chemical  reaction  (in  this  case 
bond  rupture)  is  reciprocally  dependent  on  the  exponent  of  the  energy 
barrier,  mechanisms  that  provide  a  lower  energy  activated  state  will 
enhance  crack  growth  rates.  This  brief  discussion  has  omitted  factors 
resulting  from  changes  in  the  energy  of  the  Si-O-Si  bond  due  to  adsorbed 
water  (Step  l)  and  the  energy  of  the  reaction  products  (Step  3);  these 
terms  are  not  critical  to  this  communication  and  will  be  presented  in 
detail  elsewhere. 

Experimentally,  the  role  of  the  concerted  reaction  (simultaneous 
electron  and  proton  transfer),  as  the  rate  of  limiting  step  to  fracture 
can  be  born  out  through  isotopic  effects.  Results  of  fracture  experiments 
performed  in  deuterated  vater  produced  crack  velocities  consistently 
lower  than  those  attained  in  normal  water.  This  result  can  be  analyzed 
in  terms  of  the  different  ground  state  vibrational  energies  of  deuterium 
and  prot.iura  resulting  in  different  activation  barriers  for  the  two 
isotopes.  The  lower  rack  velocities  measured  in  deuterated  vater  are 

o 

in  agreement  .t.h  uie  larger  energy  barrier  expected  for  deuteron  transfer. 

-  1.  - 


This  model  for  chemical  bond  rupture  presents  a  number  of  interesting 

implications.  First,  it  suggests  that  environments  other  than  vater 

should  enhance  cracking  in  silica  glass  if  the  species  possess  structural 

and  bonding  features  similar  to  vater  l.e.  proton  donor  sites  at  one  end 

of  the  molecule  (or  group)  and  lone  pair  orbitals  at  the  other.  The 

o 

environmental  molecule  must  also  fit  between  the  Si-0  bond  ( =1 . 63  A). 

To  test  this  hypothesis,  crack  grovth  studies  were  conducted  in  vitreous 

silica  exposed  to  ammonia  (NH^).  Ammonia  has  a  structure  similar  to  vater 

except  that  it  has  three  orbitals  o-bonded  to  hydrogens  leaving  one  lone 

pair  orbital  directed  away  from  the  hydrogen  atoms  (.the  N-H  bond  length 
o 

is  0.91*  A).  Figure  3  shows  that  indeed,  NH^  enhanced  crack  grovth  over 
the  entire  velocity  range.  When  evaluating  such  data  it  is  important  to 
be  aware  that  small  amounts  of  vater  in  an  inert  carrier  (N^)  can  have 
large  effects  on  crack  grovth  rates.  However,  as  previously  discussed, 
the  enhancement  due  to  small  amounts  of  vater  is  limited  in  velocity. 

Since  the  NH^  gas  used  in  this  experiment  contained  less  HgO  than  the 
N2  +  HgO  gas  mixture,  the  absence  of  a  plateau  indicates  that  the 
was  controlling  the  fracture  behavior.  Crack  grovth  experiments  in 
hydrazine  (N^H^ )  and  formsmide  (CH^NO).  shoved  that  these  environments 
also  controlled  crack  grovth  over  the  entire  velocity  range.  However, 
these  environments  were  less  effective  (V-K^  curve  shifted  to  higher  K^) 
then  either  vater  or  ammonia  possibly  due  to  steric  hinderence  and/or 
delocalization  of  lone  pair  orbitals.  It  is  also  noteworthy  that  other 
non-aqueous  environments  such  as  carbon  monoxide  (contains  lone  pair 
orbital  but  no  proton  donor  site)  or  acetonitrile  (has  high  dipole 
moment,  3-92D,  but  no  localized  lone  pair  orbital)  have  no  measureable 
effect  on  crack  growth  rates. 

The  model  also  suggests  that  environmentally  enhanced  crack  grovth 
of  this  type  can  only  occur  in  solids  containing  structural  bonds  between 


proton  acceptors  and  electron  acceptors.  Crack  grovth  measurements  in 
single  crystal  Si  suggest  little  if  any  effect  of  water  on  crack  grovth 
rates, ^  which  is  in  agreement  with  this  hypothesis. 

Finally,  information  suggested  hy  this  model  can  conceivably  he 
applied  to  calculations  of  reaction  rates  for  the  bond  rupturing  process. 
In  order  to  make  such  calculations  (either  by  semi-empirical  or  ab- 
initio  methods)  geometries  of  the  reactants  and  activated  complex  must 
be  known.  The  most  difficult  procedure  involves  determination  of  the 
activated  complex  since  it  does  not  exist  for  finite  times  and  cannot  be 
measured  spectroscopically.  Our  model  of  the  bond  rupturing  process 
does  suggest  a  possible  geometry  for  the  activated  complex  (Step  2  of 
Fig.  2)  which  may  be  useful  in  such,  calculations.  Further  details 
concerning  such  calculations  and  their  application  to  predictions  for 
fracture  behavior  will  be  presented  elsewhere. 
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Figure  1 
Figure  2 

Figure  3 


Effect  of  water  on  slow  crack  growth  in  vitreous  silica  (curves 
represent  best  fit  to  experimental  data). 

Schematic  representation  of  the  proposed  reaction  between  water 
and  a  strained  Si-O-Si  bond  at  the  crack  tip.  Reaction  steps 
involve;  (l)  adsorbtion  of  water  to  Si-0  bond,  (2)  concerted 
reaction  involving  simultaneous  proton  and  electron  transfer 
and  (3)  formation  of  surface  hydroxyls. 

Comparison  of  the  effects  of  water  and  ammonia  on  slow  crack 
growth  in  vitreous  silica  at  room  t emperature . 
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Effect  of  Corrosion  Processes  on  Subcrltlcal  Crack  Growth  in  Glass 
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Abstract 

Crack  growth  studies  were  conducted  on  soda  lime  silica,  soda 

borosllicate  and  two  binary  soda  silica  glasses  immersed  In  solutions  of 

1  Molar  LI*,  1  Molar  Cs*  or  deionized  water  at  different  pM  values.  A 

definite  effect  of  the  LI*  and  Cs*  was  observed  on  the  V-Kj  curves  In 

all  but  the  soda  lime  glass.  A  plateau  In  crack  velocity  In  the  range 
-8  -10 

10  to  10  m/sec  was  measured  on  the  binary  soda-silica  glasses  for 
Kj  <  0.35  MPam^*.  These  data  are  analyzed  In  terms  of  both  the  Ion 
exchange  and  SiO^  dissolution  steps  of  the  corrosion  process.  A  model 
of  crack  growth  in  corrosive  conditions  Is  proposed. 


*Now  at  the  Vitreous  State  Laboratory,  Catholic  University  of  America, 
Washington,  D.C.  20064 


Effect  of  Corrosion  Processes  on  Subcritical  Crack  Growth  In  Glass 
C.  J.  Simmons  and  S.  W.  Frelman 

Introduction 

Logarithmic  plots  of  crack  velocity  versus  stress  intensity,  ICj, 

for  glasses  in  moist  environments  typically  have  three  regions  of  behavior, 

a  portion  with  a  slope  of  10  to  40,  a  nearly  flat  portion,  and  a  steep, 

water  independent  part.  Crack  growth  in  region  I,  •vl0'*^-10“*  m/sec.  Is 

governed  by  the  rate  of  reaction  of  the  material  with  water  In  the 

environment^.  It  is  to  this  region  that  this  study  Is  addressed.  It  Is 

observed  that  the  slope  and  position  of  the  crack  velocity  vs.  Kj  curves  In 
'  o 

region  I  are  affected  by  glass  composition,  the  partial  pressure  of 

water  in  other  media,  and  the  pH  of  aqueous  solutions*.  However,  a 

more  direct  link  between  corrosion  rates  of  glasses  and  crack  growth 

behavior  has  not  been  established. 

Glass  corrosion  occurs  by  two  processes:® 

1)  Exchange  of  alkali  ions  in  the  glass  with  hydrogen  Ions 
from  the  enviornment. 

2)  Hydrolysis  of  the  Si-0  bonds  leading  to  dissolution  of 
the  silica  network. 

The  step  which  controls  the  overall  corrosion  rate  depends  on  the  glass 

composition  and  the  environment,  particularly  pH. 

The  purpose  of  this  paper  Is  to  show  that  crack  growth  behavior  In 

some  glasses,  particularly  the  binary  alkali  silicate  compositions,  can 

be  strongly  Influenced  by  alkali  ion  exchange.  The  presence  of  a  hitherto 

-8  -10 

undiscovered  crack  velocity  plateau  in  the  10  to  10  m/sec  range  for 
these  binary  glasses  will  be  shown  to  occur.  A  model  of  subcritical 
crack  growth  which  qualitatively  explains  the  position  and  slopes  of  V- 
Kj  curves  for  different  glasses  will  be  discussed. 
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Table  I 


Glass 

c n 

rv? 

Glass  Composition  (Mole  t) 
Na20  BgOg  CaO 

MgO 

tm 

CC) 

ta 

<*c) 

Soda-lime- 

* 

silica 

72 

14 

7 

4 

Soda 

boroslllcate 

70 

7 

23 

1350 

550 

33  Na.O  - 
67S102 

67 

33 

1350 

500 

25  Na20  - 

75S102  - 

75 

25 

1400 

500 
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Experimental  Procedure 


A  commercial  soda-1 ime-sllica  glass  and  three  experimental  glasses 
were  investigated.  Their  compositions  are  given  in  Table  1,  along  with 
the  melting  and  annealing  temperatures  for  the  experimental  glasses. 

The  glasses  were  prepared  from  reagent  grade  chemicals  melted  in  a 
platinum  crucible  and  stirred  during  firing.  Plates  were  cast  on  to 
steel  platens  and  annealed.  Double  cantilever  beam  specimens  50  mm  x 
25mm  x  2  mm,  were  cut  from  the  glass.  A  groove  approximately  1  m  deep 
was  along  the  center  of  *«ch  specimen  to  guide  the  crack.  The 
specimens  were  then  reaniealed.  Specimens  of  the  two  binary  soda  silica 
compositions  were  stereo  u\  a  dessicator  because  of  the  hydroscopic 
nature  of  these  glasses.  Crack  propagation  tests  were  run  using  the 
double  cantilever  beam  arrangement  as  described  by  Wiederhom.  ^  The 
specimens  were  dead-weight  loaded;  the  crack  position  was  monitored 

% 

using  a  traveling  microscope.  Stress  intensity  factors  were  calculated 
using  Equation  1: 


KI  = 


PL 


(aw) 


Yf2~yZ  (3.467  +  2.315  t/L) 


(1) 


where  P  is  the  load,  L  the  crack  length  measured  from  the  points  of  load 

application,  a  the  web  thicknesses,  w  the  total  thickness,  and  t  the 

half  width.  Tests  were  run  either  in  an  aqueous  solution  or  in  air  of 

22°C  and  ^30%  relative  humidity  (RH). 

The  primary  solutions  used  in  this  investigation  consisted  of  1 

molar  Cs+  and  1  molar  Li+  ions  in  deionized  water.  The  starting  materials 

* 

were  high  purity  UNQ^  and  CsNOj.  The  pH  of  the  solutions  was  varied 
by  the  addition  of  small  amounts  of  HC1  or  NH^OH,  and  monitored  throughout 
the  course  of  the  crack  growth  experiment. 


* 

Cerac  Co. 


Results 


Figure  1  shows  crack  velocity  as  a  function  of  Kj  for  soda-lime 

glass  in  the  1  M  11*  and  1  H  Cs*  solutions  at  three  different  pH  values. 

Within  the  scatter  of  the  data,  one  cannot  differentiate  between  the 

crack  growth  behavior  in  the  different  solutions.  This  is  not  entirely 

unexpected  since  Wiederhom  and  Johnson*  found  effectively  no  change  in 

crack  growth  behavior  with  changes  in  buffered  solutions  of  pH  1.8  to  9 

for  soda-lime  glass.  One  point  to  note  is  the  apparent  crack  arrest  at 

a  Kj  of  -v.27MPam1/2  for  the  soda-lime  glass  in  deionized  water.  When 

the  apparently  stationary  crack  was  reloaded  to  higher  Kj  it  began  to 

propagate  at  the  velocity  expected  from  the  crack  growth  curves.  The 

change  to  a  steeper  slope  of  the  pH  3.5  V-Kj  curves  at  lower  velocity  Is 

*  2 

in  agreement  with  the  data  reported  by  Wiederhorn  and  Bolz  for  crack 
growth  in  water  at  various  temperatures.  While  there  is  the  suggestion 
of  a  change  to  a  lower  slope  in  the  pH  9  curve,  there  is  not  enough  data 
to  be  definitive. 

Borosilicate  glass  (Figure  2)  appears  to  behave  differently  when 
tested  in  Cs+  and  Li+  solutions.  For  any  given  pH  the  data  taken  In  Cs+ 
solution  has  the  same  slope  but  lies  at  higher  Kj  than  does  the  data  In 
Li+.  The  pH  7  and  3.5  data  overlap  while  that  for  pH  9  lies  at  higher 
stress  intensities.  One  interesting  feature  of  crack  growth  in  this 
glass  is  the  crack  arrest  at  Kj  ^  0.35  MPam^2.  The  behavior  of  the 
cracks  upon  reloading  differed  markedly  from  that  In  the  soda-lime 
glass;  loading  the  arrested  cracks  to  higher  Kj  values  did  not  cause 
them  to  repropagate.  As  the  cracks  were  loaded  to  higher  and  higher  Kj 
values  catastrophic  failure  finally  occurred.  Similar  behavior  was 

4 

observed  by  Wiederhom  and  Johnson  in  Pyrex  tested  In  acid  environments. 
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In  this  borosllicate  glass  it  is  impossible  to  avoid  phase-separation 
during  cooling  of  bulk  samples;  additional  microstructure  growth  occurs 
during  the  anneal  at  550°C.  The  samples  tested  had  an  interconnected 
microstructure  (5  -  10  NM)  of  equal  volume  fraction  of  which  one  phase 
contains  <50%  Si Og .  This  phase  is  less  resistant  to  chemical  attack 
than  any  of  the  other  glasses  tested  and  may  account  for  the  severe 
crack  blunting  behavior  observed  only  in  this  glass.  Crack  blunting  is 
further  indicated  by  the  appearance  of  ridge  marks  emanating  from  the 
position  of  the  stopped  crack. * 

Crack  growth  data  for  the  33%  soda-silica  glass  is  presented  in 
Fig.  3.  Because  of  the  many  curves,  the  individual  data  points  are  not 
shown.  There  are  several  important  features  in  this  plot.  A -marked 
difference  between  the  V-Kj  curve  measured  in  air  of  30%  RH  is  its  location 
at  significantly  higher  Kj  and  greater  slope  compared  to  those  measured  in 

+  4. 

aqueous  solutions.  The  V-Kj  data  in  Cs  and  Li  solutions  is  clearly 

separated.  In  contrast  to  the  borosilicate  glass,  the  data  taken  in  IN 

li+  lies  at  higher  Kj  values  than  does  the  Cs+  data  at  the  same  pH,  and 

also  has  a  significantly  higher  slope.  The  effects  of  pH  overlap  this 

behavior,  with  the  data  for  pH  6.5  lying  at  slightly  lower  Kj  values 

than  either  pH  3.5  or  pH  9.  One  significant  point  is  that  the  data  in  IN 

Cs+  is  identical  with  that  for  deionized  water.  The  plateau  at  low  Kj  observed 

in  this  glass  can  be  seen  to  be  unaffected  by  the  presence  of 

either  Cs+  or  Li+,  but  the  crack  velocity  at  which  the  plateau  growth 

occurs  is  dependent  on  the  pH  of  the  solutions,  with  pH  3.5  lying  at 

higher  velocities  than  pH  6.5.  This  plateau  has  been  found  to  extend  to 

1/2 

Kj's  as  low  as  0.04  MPam  .  Data  taken  at  pH  9  does  not  show  a  distinct 
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plateau;  crack  velocity  continues  to  decrease  with  decreasing  ,  but 
decreases  at  a  much  lower  rate  than  in  the  other  portion  of  the  curve. 

The  crack  growth  behavior  in  the  25*  soda-silica  glass*  Figure  4* 

is  quite  similar  to  that  in  the  other  binary  glass,  namely^the  data  in 

30X  RH  air  lies  at  higher  Kj  than  does  the  data  in  the  Li*  and  Cs* 

solutions.  One  interesting  feature  of  the  data  for  this  glass  is  the 

definite  indication  of  crack  arrest  in  air  in  direct  contrast  to  the 

oppositely  shaped  curves  in  the  Li+  and  Cs+  solutions.  The  difference 

4-3 

in  crack  velocity  level  in  solutions  with  and  without  200  ppm  of  A1  is 
related  to  the  dissolution  characteristics  of  this  glass  as  discussed  in 
the  next  section. 

Discussion 

It  is  suggested  that  the  crack  growth  behavior  described  previously 
can  be  explained  by  a  balance  between  stress  enhanced  chemical  reactions 
and  the  normal  corrosion  of  these  glasses  in  the  different  solutions.  * 
The  overall  corrosion  rate  of  the  glasses  studied  increases  in  the 
order,  soda-lime-silica,  borosllicate,  25%  soda-silica,  33*  soda-silica. 
The  soda-lime-silica  glass  is  relatively  durable,  so  a  crack  growth 
behavior  essentially  independent  of  the  exchangeable  species  in  solution, 
i.e.  Li+  and  Cs,  and  pH  between  3.5  to  9.0  is  not  unexpected. 

In  the  borosilicate  glass,  the  effect  of  the  particular  ions  in 
solution  is  reasonable  clear.  As  noted  in  the  introduction,  one  step  in 
the  glass  corrosion  process  is  the  exchange  of  ions  in  the  glass  for 
ions  in  solutions,  i.e.,  normally  H+.  In  this  case  some  exchange  of  Li+ 
and  Cs+  ions  for  Na+  apparently  occurs. 
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In  very  non-durable  glasses  such  as  the  binary  soda-silicates, 
drastic  effects  of  changing  the  corrosive  solution  might  be  expected. 

In  this  case  Cs+  seems  to  produce  the  same  crack  growth  behavior  as 
water.  It  Is  likely  that  Cs+-Na+  and  H+-Na+  exchange  rates  are  comparable 
at  these  solution  concentrations  and  that  glasses  containing  large 
concentrations  of  either  Cs+  or  H+  exhibit  similar  corrosion  behavior. 
Evidence  for  the  latter  Is  the  similar  corrosion  rates  measured  by 
infrared  reflection  spectroscopy  for  NaOg-SIOg  glasses  In  either  1M  Cs+ 
solutions  or  deionized  water.®  In  L1+  solutions  the  opposite  is  true; 
lithia-silica  glasses  are  more  durable  than  soda-silica  glass,®  so  that 
any  exchange  of  Li+  for  Na+  tends  to  produce  a  protective  layer  which 
would  retard  the  rate  of  crack  growth.  It  is  hypothesized  that  the 
small  changes  seen  in  the  curves  with  variations  in  pH  may  reflect 
variations  in  the  degree  of  hydration  of  Cs+  and  J.1+  and  therefore 

4.  ' 

reflect  their  ability  to  exchange  with  the  Na  ions  in  the  glass.  One 
possibility  which  cannot  be  ignored  but  which  will  not  be  dealt  with  in 
this  paper  is  that  ion  exchange  of  H+  for  Na+  produces  stresses  at  the 
crack  tip  which  influence  the  growth  behavior.  Evidence  for  this  is  the 
higher  plateau  velocity  at  lower  pH  plus  the  observation  of  small  cracks 
in  the  bulk  glass  at  pH-1.  It  was  also  observed  that  some  peeling  of 
specimen  surfaces  occurred  after  drying. 

When  Kj  is  lowered  to  the  point  where  the  rate  of  stress-aided- 
dissolution  equals  the  normal  rate  of  corrosion,  a  change  in  crack  tip 
radius  would  be  expected, as  predicted  from  the  theory  of  Charles  and 
Hillig.^  In  the  soda-1 ime-sil lea  and  especially  the  borosilicate  glass, 
this  change  in  crack  tip  radius  leads  to  crack  arrest.  A  different 
phenomenon  occurs  in  the  soda  silica  glasses,  namely  the  crack  continues 
to  grow  at  a  velocity  which  approximates  the  rate  of  silica  dissolution 
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in  these  glasses.  The  latter  hypothesis  is  based  on  rates  of  silica 

dissolution  calculated  from  the  work  of  Douglas  and  El  Shany®  who  measured 

silica  extraction  rates  which  would  correspond  to  removal  of  silica 

-9  -8 

layers  in  the  range  of  10  -10  'm/sec.  The  saturation  of  the  solution 
at  the  crack  tip  with  silica  at  pH  9  may  explain  the  absence  of  an 
observable  plateau  under  these  conditions.  However,  a  question  arises 
as  to  why  these  cracks  continue  to  grow  rather  than  round  out  and  stop. 

It  is  speculated  that  growth  continues  because  the  surface  that  forms, 
either  as  a  result  of  dealkilization  and/or  reprecipitation  of  dissolved 
material,  is  mechanically  weak  and  relatively  nonadherent.  It  Is  proposed 
that  even  a  small  stress  at  the  crack  tip  produces  a  break  In  the  corrosion 
layer  allowing  fresh  solution  to  reach  the  unreacted  glass.  The  small 
stress  aids  removal  of  material  at  the  tip  in  preference  to  removal  at 
other  points  on  the  crack  surface,  so  that  comple'te  blunting  does  not, 
occur. 

If  crack  growth  rates  in  this  plateau  are  governed  by  the  rate  of 
stress  assisted  silica  dissolution,  one  would  expect  that  changes  In  the 
solutions  which  lead  to  changes  in  corrosion  rates,  would  be  mirrored  in 

o 

changes  in  this  crack  growth  plateau.  Dilmore,  Clark  and  Hench  showed 

-f-3 

that  the  addition  of  200  ppm  of  A1  ions  to  the  solution  lowers  the  silica 
dissolution  rate  in  a  binary  lithia-silica  glass  by  about  an  order  of 
magnitude.  As  shown  in  Fig.  4  for  the  252  soda  silica  glass,  the  addition 
of  200  ppm  Al+3  also  lowers  the  crack  velocity  plateau  In  distilled  H^O 
by  an  order  of  magnitude,  giving  further  credence  to  the  idea  of  this 
plateau  being  strongly  influenced  by  the  S102  dissolution  rate.  One 
point  which  is  unclear  is  why  the  252  and  332  Na20  glasses  have  plateaus 
in  the  same  range,  when  their  $102  dissolution  rates  should  be  quite 
different.  One  possibility  is  that  the  solution  at  the  crack  tip  becomes 
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saturated  with  Si 02  and,  therefore,  stress  assisted  dissolution  cannot 
occur  faster  than  the  SlOg  diffusion  rate  away  from  the  crack  tip. 

This  concept  of  a  change  from  stress-corrosion  controlled  crack 
growth  to  one  involving  mainly  dissolution  implies  that  the  crack  shape 
should  also  change.  If  the  crack  shape  changes  in  going  from  region  I 
to  the  low  level  plateau  or  vice-versa,  one  should  be  able  to  detect  the 
rounding  or  sharpening  of  the  crack  with  time.  In  order  to  check  this, 
an  experiment  was  conducted  in  which  a  crack  in  the  25*  soda  silica 

glass  immersed  in  a  1M  Li  solution  with  200  ppm  A1  (pH*6.5)  was 

-9  1/2 

allowed  to  reach  a  steady  state  velocity  of  *10  m/sec.  at  a  Kj  of  0.15  HPam  '  . 

The  stress  intensity  was  abruptly  raised  to  0.50MPam^Z,  at  which  point 

-5 

the  steady  state  crack  velocity  should  have  been  *10  m/sec.  Even  though 
an  initial  sharp  increase  in  V  occurred,  (Fig.  5),  the  crack  velocity 
reached  the  expected  value  only  after  approximately  1  hour  under  load,  after 
which  time  the  velocity  remained  constant.  The  measurement  of  the 
initial  jump  could  be  in  considerable  error  because  of  the  small  amount 
of  crack  extension  that  occurred.  This  data  suggests  that  a  gradual 
reduction  in  crack  tip  radius  was  required  (the  balance  between  corrosion 
and  stress  corrosion  had  to  be  readjusted)  before  a  steady  state  velocity 
could  be  reached  at  the  higher  Kj. 

Variations  in  crack  tip  shape  are  postulated  to  occur  throughout 
the  regime  of  subcritical  crack  growth.  From  the  data  for  region  I 
crack  growth  (Fig.  6)  one  can  see  a  significant  difference  between  the 
slopes  and  positions  of  the  different  curves  depending  on  the  glass 
composition  and  test  environment.  If  one  looks  in  particular  at  the 
binary  soda-silica  glasses,  it  is  apparent  that  the  V-K^  curves  for  the 
less  chemically  durable  33*  soda  glass  lie  at  higher  Kj  values  than  the 
25*  soda  glass  in  either  Cs+  or  L1+solutions.  This  is  opposite  to  what 
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one  might  initially  expect  from  the  rates  of  corrosion.  This  feature  of 

the  data  suggests  that  crack  tip  shapes  may  differ  in  different  glass 

compositions  and  that  the  less  durable  glasses  have  larger  crack  tip 

radii  and  therefore  smaller  stress  fields  {than  would  be  calculated  from 

fracture  mechanics  theory)  at  their  tips.  One  possible  model  would 

suggest  that  the  two  primary  factors  determining  the  position  and  slope 

of  a  V-Kj  curve  are  the  rate  of  corrosion,  a  process  which  tends  to 

round  the  crack  tip,  reducing  the  stress  intensity,  and  the  stress 

dependence  of  the  corrosion  process  which  tends  to  produce  a  sharper 

crack.  This  concept  is  in  direct  analogy  to  the  model  of  Charles  and 

Hi 11 ig. ^  We  are  suggesting  that  each  glass  composition  has  its  own 

crack  tip  radius  which  decreases  with  increasing  stress  intensity  and 

that  this  decrease  in  crack  tip  radius  with* Increasing  stress  Intensity 

gives  rise  to  the  slope  in  the  V-Kj  curves.  One  anomally  to  this  model 

is  fused  silica,  which  should  have  a  very  low  corrosion  rate,  yet  He* 

well  to  the  right  of  any  of  the  crack  growth  curves  shown  in  Figure  6. 

It  may  be,  however,  that  the  type  of  corrosion  product  formed  may  also 

have  a  significant  influence  on  the  crack  growth  behavior  so  that  one 

may  need  to  compare  families  of  glasses  rather  than  compare  glasses  of 

widely  varying  composition.  This  kind  of  comparison  Is  suggested  by 

g 

the  work  of  Hench  et.  al.  who  showed  that  corrosion  behavior  in  glasses 

could  be  of  at  least  5  different  types.  Another  possible  anomally  is 

the  behavior  of  the  soda-silica  glasses  in  air,  where  their  resistance 

to  stress  corrosion  seems  to  be  significantly  higher,  even  though  one  would  expect 

the  crack  tip  radii  to  be  smaller  because  of  lower  corrosion  rates. 

However,  corrosion  products  may  be  easily  removed  in  solutions  while  being 
trapped  at  the  crack  tip  in  air,  thereby  hindering  the  reactants  from 
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reaching  the  unreacted  material.  Secondly,  there  is  a  much  smaller 
concentration  of  water  to  react  with  the  surface.  Both  factors  would 
lead  to  reduced  stress  corrosion  rates. 


Summary  and  Conclusions 

Crack  velocities  have  been  measured  as  a  function  of  Kj  for  four 
different  glasses  in  solutions  of  Li+  and  Cs+  held  at  different  pH.  The 
data  for  soda  borosilicate  and  the  two  binary  soda-silica,  glasses  show 
definite  effects  of  the  presence  of  Li+  or  Cs+  ions.  In  the  binary 
glasses,  the  V-Kj  curves  in  L1+  lie  at  higher  Kj's  because  of  the  formation 
of  a  protective  layer,  and  are  steeper  than  those  obtained  in  either  Cs+ 
or  deionized  water  because  the  ion  exchange  rate  of  Li+  for  Na+  is  slow 
enough  so  that  as  the  crack  velocity  increases,  the  concentration  of  L1+ 
in  the  surface  is  steadily  reduced  until  the  protective  layer  no  longer 
forms.  In  the  case  of  Cs+  and  H+,  ion  exchange  rates  are  higher  and  so 
can  keep  up  with  crack  growth  rates  over  a  much  larger  range  of  velocities. 
Changes  in  pH  are  shown  to  produce  smaller  shifts  in  these  curves. 

A  plateau  in  crack  velocity  is  observed  for  the  binary  glasses  at 
0.04  <  Kj  <  0.35  MPam1^2.  The  value  of  the  velocity  plateau  appears  to 
depend  on  the  SiC^  dissolution  rates  of  the  glass  in  the  given  solution, 
but  some  affect  of  local  stresses  due  ion  exchange  cannot  be  ignored. 

It  is  concluded  from  crack  acceleration  experiments  that  changes  in 
crack  tip  shape  must  take  place  as  a  crack  propagates  from  the  low  level 
plateau  into  Region  I.  A  model  of  environmentally  assisted  crack  growth 
is  put  forth  which  suggests  that  the  position  and  slope  of  the  V-Kj 
curve  for  a  particular  glass  in  a  given  environment  is  governed  by  a 
balance  between  the  tendency  of  corrosion  processes  to  increase  the 
crack  tip  radius  and  the  stress  corrosion  process  which  would  sharpen 
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Figure  Captions 


Figure  1 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Crack  Velocity  -  Kg  plots  for  soda  lime  silica  glass  tested  In 
IN  Li*  and  1M  Cs*  solutions,  and  deionized  water.  Note  apparent 
cracks  arrest  at  Kj  ■  0.27  MPam^  slope  of  the  line  through  the 
data  *  19. 

Crack  velocity  -  Kj  plots  for  soda  borosilicate  glass  tested  in 
the  same  solutions  as  in  Figure  1.  Note  the  definite  separation 
in  the  curves  for  Cs*  and  Li*,  with  the  former  lying  at  higher  Kg's  than 
the  latter  for  all  pH  ranges  investigated.  The  slope  of  all 
the  curves  is  *  35. 

Crack  velocity-Kj  plots  for  33%  Na20-S1Ck  glass  In  the  different 
environments.  Note  that  curves  -in  LI*  are  steeper  and  higher  Kj 
than  those  in  Cs*.  The  data  taken  in  deionized  water  overlaps  that 
in  Cs*. 

Crack  velocity  -Kj  plots  for  25%  Na20-S102  glass  In  different 
environments.  Note  that  data  obtained  In  any  solution  containing 
~200  ppm  At  lies  on  the  lower  of  the  two  plateaus.  A  definite 
crack  arrest  was  observed  in  this  glass  tested  In  air. 

Crack  velocity  In  25%  Na20-S102  glass  as  a  function  of  time  at  a 

stress  intensity  of  0.53  MPam^.  Crack  was  initially  propagating 

.0 

at  a  velocity  of  ~10  m/sec.  when  Kj  was  raised. 

Comparison  of  V-Kj  curves  in  Region  I  for  all  glasses 
Investigated  In  this  study. 
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INTENSITY 


Raman  investigation  of  optical  fibers  under  high  tensile  stress  a) 

G.  E.  Walrafen 

Department  of  Chemistry,  Howard  University,  Washington,  D.  C.  20059 

P.  N.  Krishnan 

Department  of  Chemistry,  Coppin  State  College,  Baltimore,  Maryland  21216 

S.  W.  Freiman 

National  Bureau  of  Standards,  Gaithersburg,  Maryland  20234 
(Received  29  July  1980;  accepted  for  publication  10  November  1980) 

Raman  spectra  have  been  obtained  from  fused  silica  optical  fibers  under  tensile  stresses  from  0  to 
3.3  GPa  (33  kbars).  Reversible  intensity  increases,  relative  to  the  principal  Raman  maximum  at 
440  cm  ~ 1 ,  were  observed  for  the  defect  peak  at  490  cm  “ 1  and  for  shoulders  near  350-375  and  1 1 5 
cm  " 1 .  Application  of  tensile  stress  to  fused  silica  appears  to  produce  changes  in  the  stretched  Si — 
O  defect  sites,  as  well  as  changes  in  the  main  network  structure. 

PACS  numbers:  63.50.  +  x,  78.30.Jw,  81.20.Pe 


INTRODUCTION 

Raman  features  from  fused  silica  near  490  and  604 
cm  " 1  have  been  examined  in  several  recent  investigations, 1-7 
see  Fig.  1  (stars).  Raman  intensities  at  490  and  604  cm" ' 
were  found  to  increase  with  increasing  Active  temperature 


Ty  when  the  OH  content  was  constant  (or  zero).4  The  490- 
and  604 -cm  ~ 1  Raman  intensities  were  also  observed  to  de¬ 
crease  with  increasing  OH  content  at  constant  Ty.*  In  the 
latter  case  a  weak  Raman  band  near  970  cm" due  to  Si 
versus  OH  stretching,4  "  was  produced  at  the  expense  of  in- 


800  600  400  200 


FIG.  I.  Forward  Raman  scattering  from 
fused  silica  optical  fiber  stressed  to  3.34 
GPa  |  upper),  and  unstressed  (lower).  Spec¬ 
tra  in  the  inset  were  obtained  in  the  se¬ 
quence  unstressed  U ;  stressed  5  to  2. 19 
GPa;  and  relaxed  R.  Note  that  the  inset 
spectra  have  been  moved  to  the  left — their 
position  does  not  correspond  to  the  cm  "  ' 
values  below  them.  Baselines  A  and  B  are 
shown  by  dashes. 
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tensity  at  the  490-  and  604-cm~ 1  positions.  Other  findings, 
such  as  the  observation  that  the  490-  and  604-cm" 1  intensi¬ 
ties  from  samples  having  high  TF's  could  be  reduced  by  an¬ 
nealing,  strengthened  the  hypothesis  that  these  peaks  could 
be  assigned  to  defects  in  the  silica  structure.4 

The  nature  of  the  above  defects  suggested  that  their 
concentration  and/or  molar  intensity  should  increase  under 
the  application  of  tensile  stress.  In  order  to  subject  a  suffi¬ 
cient  volume  of  glass  to  large  tensile  stresses,9  experiments 
were  conducted  on  glass  optical  fibers.  An  interpretation  of 
the  changes  observed  in  the  Raman  spectrum  is  presented  in 
terms  of  variations  in  the  Si — O  binding  at  defect  sites. 

EXPERIMENTAL  PROCEDURES 

Pure  fused  silica  optical  fibers  were  used  in  all  experi¬ 
ments  (silicone  rubber  cladding,  200  /im  diameter).10  The 
OH  content  of  the  fibers  was  about  530  ppm, 1 1  and  TF  values 
were  estimated  to  be  roughly  1700-1800  *C. 12  The  constan¬ 
cy  of  the  OH  concentration  was  assured  by  monitoring  the 
Raman  OH-stretching  intensity  near  3700  cm”1.13 

Forward  scattered  Raman  spectra  were  obtained  by  fo¬ 
cusing  the  fiber  output  on  the  entrance  slit  of  an  Instruments 
S.  A.  HG2S  holographic  grating  double  monochromator 
(Fig.  2).  Raman  excitation  was  accomplished  using  5 14.5-nm 
argon  ion  laser  radiation  at  about  0.6  W.  Slit  widths  ranged 
from  2  to  3  cm- ',  and  detection  was  accomplished  with  an 
uncooled  Hamamatsu  R928  PM  tube,  a  Keithley  414S  pi- 
coammeter,  and  a  high-speed  Esterline  Angus  LI  101 S  re¬ 
corder. 

Load  was  applied  to  the  fiber  by  suspending  a  weight 
from  a  clamp  near  the  end  of  the  fiber  (Fig.  2).  The  optical 
fiber  was  held  by  strips  of  sponge  rubber  placed  between 
wooden  strips  30-50  cm  in  length,  tightly  clamped  together 
in  several  places.  This  arrangement  prevented  any  slippage 
during  application  of  load.  The  stressed-to-unstressed  length 
ratio  was  maintained  in  the  range  of  5-50  to  reduce  the  un¬ 
stressed  spectral  contribution. 

Some  Raman  experiments  were  also  attempted  in 
which  only  stressed  fiber  was  examined  by  focusing  90°  scat¬ 
tered  radiation  on  the  slit,  but  these  experiments  were  aban¬ 
doned  because  the  S/N  (signal-to-noise)  ratio,  although 
good,  was  at  least  20-fold  lower  than  that  obtained  from 
forward  Raman  scattering. 

Because  the  Raman  depolarization  ratio  varies  across 
the  fused  silica  spectrum,14  it  was  necessary  to  determine  if 
the  polarization,  which  is  ordinarily  scrambled  by  an  optical 
fiber,  changes  under  tensile  stress.  Preliminary  experiments 
involved  examination  of  the  mode  structure  and  the  polar¬ 
ization  at  the  fiber  output.  No  significant  changes  were  ob¬ 
served  in  the  high-order  mode  structure  under  stress,  and  no 
changes  whatever  in  polarization  were  detected  using  polar¬ 
izers  or  crossed  polarizers. 

Quantitative  measurements  of  the  polarization  were 
made  using  the  optical  fiber  emission  to  excite  both  forward 
and  90°  Raman  scattering  in  liquid  CCl4.  Ratios  of  integrat¬ 
ed  Raman  intensities  were  determined  for  the  CC14  peaks  at 
218(dp),  314(dp),  and  459(p)  cm  '.  (Here  dp  refers  to  a  de¬ 
polarization  ratio  of0.75and  ptoaratio  of  about  0.01. I*  For 


FIG.  2.  Schematic  illustration  of  method  used  to  obtain  forward  Raman 
spectra  (him  an  optical  fiber  while  applying  tensile  stress. 

both  forward  and  90°  scattering,  no  changes  whatever  were 
detected  in  the  intensity  ratios  /45<//3i4  and  /45»//21,  for  ten¬ 
sile  stresses  to  2.19  GPa.  Further,  the  intensity  ratios  ob¬ 
tained  from  CCI4  were  those  characteristic  of  unpolarized 
excitation,  as  expected  from  the  fact  that  polarization  is 
scrambled  by  internal  reflection  in  high-mode  optical  fibers. 
Also,  it  should  be  emphasized  that  the  above  integrated  Ra¬ 
man  intensity  ratios  from  CC14  constitute  a  particularly  sen¬ 
sitive  test  of  the  polarization  properties  of  the  exciting  radi¬ 
ation  emitted  from  the  optical  fiber. 

RESULTS 

Stressed  and  unstressed  Raman  spectra  illustrative  of 
the  approximately  50  individual  experiments  performed  are 
shown  in  Fig.  1.  The  (upper)  spectrum  was  taken  at  a  tensile 
stress  of  3.34  GPa  (33.4  kbars),  while  the  lower  spectrum  was 
taken  from  an  unstressed  fiber.  The  intensification  of  the 
490-cm_  1  defect  peak  (starred),  relative  to  the  440-cm" 1 
peak,  is  visually  obvious,  and  is  emphasized  by  the  horizon¬ 
tal  lines.  In  the  unstressed  spectrum,  the  490-  and  440-cm  “ 1 
peaks  occur  at  the  same  vertical  height  as  seen  by  the  hori¬ 
zontal  line  through  them,  but  at  3.34  GPa  (upper)  the  490- 
cm  ' 1  peak  height  is  well  above  the  horizontal  line.  Similar 
conclusions  may  be  drawn  from  examination  of  the  inset  of 
Fig.  1  which  refers  to  three  spectra  run  in  the  sequence  un¬ 
stressed  U;  stressed  S  at  2. 19  GPa,  S  and  relaxed  R  with  20- 
30  min  between  spectra.16 

The  visual  observation  that  the  490-cm  ~ 1  peak  height 
increases  relative  to  the  440-cm  ~ 1  peak  height,  i.e.,  that 
/4<*//44«  increases,  may  mean  either  that  /49ll  increased  at 
constant  /44„,  or  that  /440  decreased  at  constant  /4W,  or  that 
both  /4„,  and  /440  changed,  or  even  that  more  complicated 
changes  occurred,  such  as  broadening,  etc.  Hence,  if  /44„  is  to 
be  regarded  as  a  reference  standard,  it  is  desirable  to  under¬ 
stand  how  /440  varies  relative  to  other  peaks,  such  as  the 
intense  peak  at  60  cm  ~ 

Peak  heights  above  baseline  A  in  Fig.  1,  were  measured 
for  the  60-,  440-  and  800-cm  ~ 1  peaks  at  a  series  of  tensile 
stresses.  Consistency  of  drawing  baseline  A  under  all  spectra 
was  obtained  by  using  a  single  template  whose  shape,  base¬ 
line  A  in  Fig.  1 ,  was  determined  through  repeated  trials.  The 


-  34  - 


2633 


J.  Appl.  Phys.,  Vol.  52.  No.  4,  April  1981 


Walrafen,  Krishnan,  and  Freiman 


2833 


template  was  positioned  to  a  given  Raman  spectrum  by  si¬ 
multaneously  matching  two  regions:  (1)  the  region  between 
the  exciting  line  and  the  sharp  minimum  near  20  to  25  cm  ~ 1 
and  (2)  the  region  from  about  900-1000  cm-1.  The  results 
are  shown  in  Table  I  in  terms  of  the  ratios  /so o/h*o> 

and  /sooZ/fco- 

The  average  value  for  the  peak  height  ratio  Ah/Am 0 
from  Table  I  is  0.70  ±  0.02,  which  within  present  errors  is 
constant.  Similarly,  the  ratios  Iioo/Iuo  and  /wx/Aso  are  also 
seen  to  be  constant  (Table  I).  The  constancy  of  the  ratio 
Ax/Amo  is  particularly  important  because  the  peaks  at  60  and 
440  cm  - 1  have  both  the  largest  peak  heights  and  the  highest 
integrated  intensities  (Fig.  1)  in  the  fused  silica  spectrum. 

Quantitative  values  of  peak  height  ratios  are  often  diffi¬ 
cult  to  obtain  accurately  when  weak  features  such  as  at  490 
and  604  cm  - 1  are  involved  because  such  ratios  suffer  from 
errors  in  baseline  estimates  and  peak  S/N  ratios.  Peak  height 
ratios  also  ignore  component  broadening.  A  better  method  is 
to  use  ratios  of  integrated  intensities.  For  the  490-  and  604- 
cm-  1  peaks  (Fig.  1),  integrated  intensities  may  be  obtained 
using  baseline  B.  Baseline  B  was  obtained  by  using  a  French 
curve.  In  contrast,  an  integrated  intensity  for  the  440-cm-  1 
component  could  not  be  readily  obtained  because  the  com¬ 
ponent  shape  is  unknown.  Therefore  the  total  integrated  Ra¬ 
man  intensity  from  about  20  to  900  cm  ~ 1  was  obtained  using 
baseline  A  minus,  of  course,  the  contributions  from  the  490- 
and  604-cm  - 1  components.  The  resulting  ratios  of  integrat¬ 
ed  Raman  intensities  I490/Ir  and  1(^/1  T,  where  T  refers  to 
the  total  Raman  area  minus  the  490-  and  604-cm  ~ 1  compo¬ 
nent  areas,  are  shown  for  tensile  stresses  from  0  to  2. 19  GPa 
in  Table  II.  In  addition,  the  quantities  /49O//*40  and  Ium/I«0 
are  listed  in  T able  II,  where  f490  and  are  integrated  inten¬ 

sities  above  baseline  B,  and  /4 40  refers  to  the  peak  height 
above  baseline  A.  The  latter  ratios  were  used  in  Table  II 
because  of  the  intensification  of  shoulders  described  subse¬ 
quently.  It  should  also  be  made  very  clear  that  the  490-cm  - 1 
component  area  (the  integrated  component  intensity)  ob¬ 
tained  above  baseline  B  is  independent  of  effects  due  to 
broadening  of  the  440-cm- 1  component.  Slight  broadening 
of  the  440-cm-  ‘  component  was  dearly  evident  from 
changes  in  the  shape  of  baseline  B  under  the  490-cm  - 1  peak, 
but  this  broadening  effect  was  completely  negated  by  use  of 
baseline  B. 

It  is  apparent  from  Table  II  that  the  only  significant 
changes  involve  the  490-cm  - 1  component,  as  previously 


TABLE  1.  Raman  peak  height  ratios  (I  refers  to  peak  height)  for  the  60-, 
440-,  and  800-cm" '  peaks  from  fused  silica  for  tensile  stresses  to  2. 19  GPa. 


Siren 

(GPa) 

^0^440 

AmJo/^440 

1*0 t/l*0 

0 

0.68 

0.17 

0.25 

0.31 

0.69 

0.17 

0.24 

0.62 

0.69 

0.17 

0.25 

0.94 

0.71 

0.17 

0.24 

1.25 

0.72 

0.17 

0.24 

1.56 

0.71 

0.18 

0.25 

1*7 

0.67 

0.17 

0.26 

2.19 

0.71 

0  18 

0.25 
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TABLE  II.  Ratios  of  integrated  Raman  intensities  and  peak  heights  for  the 
440- ,  490- ,  and  604-cm  ~ '  peaks  from  fused  silica.  /4W  and  refer  to 
integrated  component  intensities  above  baseline  B  (Fig  ij.  /,  is  the  total 
integrated  intensity  above  baseline  A  minus  the  490-  and  604-cm  " 1  compo¬ 
nent  areas  above  baseline  B.  lw  is  the  peak  height  above  baseline  .4  at  440 
cm" 


Stress 

(GPa) 

/4W/ /  r 

f 490^^440 

J*0i/^44Q 

0 

0.011 

0.013 

0.0009 

0.0011 

0.31 

0.011 

0.015 

0.0009 

0.0012 

0.62 

0.012 

0.014 

0.0010 

0.0011 

0.94 

0.014 

0.013 

0.0011 

0.0011 

1.25 

0.016 

0.013 

0.0013 

0.0011 

1.56 

0.016 

0.015 

0.0013 

0.0012 

1.87 

0.016 

0.013 

0.0013 

0.0011 

2.19 

0.017 

0.013 

0.0015 

0.0011 

suggested  from  examination  of  Fig.  1 .  The  original  data  lead¬ 
ing  to  Table  II  (using  three  figures,  instead  of  the  two  signifi¬ 
cant  figures  of  Table  II)  were  treated  by  linear  least  squares. 
Changes  of  about  62%  in  /490//r.  and  of  about  67%  in 

resulted  for  stresses  to  2.19  GPa,  using  the  least- 
squares  estimates,  but  changes  in  the  ratios  involving  the 
604-cm  ~ '  component  were  not  significant  (negligible).  Be¬ 
cause  the  relative  604-cm  - 1  intensity  did  not  change  with 
tensile  stress,  it  can  be  concluded  from  Tables  I  and  II  that 
Raman  intensities  at  60,  440, 604,  and  800  cm- 1  are  either 
all  independent  of,  or  vary  in  the  same  way  with,  tensile 
stress. 

In  addition  to  relative  intensification  of  the  490-cm  - 1 
peak,  other  small  spectral  changes  were  seen  for  Raman  fea¬ 
tures  below  440  cm-1.  One  effect  which  occurred  consis¬ 
tently  throughout  this  work  is  an  apparent  enhancement  un¬ 
der  tensile  stress  of  a  shoulder  near  350-375  cm- ',  evident 
from  careful  examination  of  Fig.  I,  particularly  the  inset.  In 
the  inset,  the  region  from  about  440  to  350  cm- 1  has  the 
same  slope  in  both  the  U  and  R  spectra.  However,  this  slope 
is  smaller  for  the  S  spectrum,  indicating  a  slight  intensifica¬ 
tion  of  the  350-375  cm- 1  shoulder  component.  However, 
the  350-375  cm  - 1  shoulder  would  also  appear  to  intensify  if 
the  440-cm  - 1  component  broadened.  Experience  with  base¬ 
line  B  at  high  tensile  stress  indicates  a  slight  but  definite 
broadening  of  the  440-cm  - 1  component,  which  could  ac¬ 
count  for  part  (or  all)  of  the  intensification  at  350-375  cm-1. 
In  regard  to  intensification  due  to  broadening  of  neighboring 
components,  it  should  be  emphasized  again  that  the  relative 
intensification  at  490  cm  - 1  cannot  be  so  explained  because 
the  440-cm  - 1  broadening  is  subtracted  from  the  integrated 
intensity  at  490  cm  - 1  by  virtue  of  curvature  in  baseline  B 
(Table  II). 

Some  shape  changes  may  also  occur  in  the  Raman  re¬ 
gion  between  0  and  200  cm  ~ 1  under  increasing  tensile  stress. 
A  shoulder  at  roughly  1 15  ±  5  cm- 1  appears  to  be  more 
prominent  under  high  stress,  compared  to  the  nominal  60- 
cm  - 1  intensity  maximum  (Fig.  1).  The  Raman  amplitude 
ratio  /, using  the  methods  baseline,  increases  roughly 
5%  in  the  figure.  The  position  of  the  Raman  peak  at  60  cm- 1 
seemed  not  to  rise  in  position  significantly,  however.  For 
tensile  stresses  from  0  to  2. 19  GPa,  the  position  ranged  from 
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60  i  2  to  63  ±  2  cm-’,  i.e.,  about  constant.  However,  it  was 
not  possible  to  determine  whether  or  not  the  relative  intensi¬ 
fication  at  1 15  cm  ~ 1  resulted  from  an  increasing  half-width 
of  the  60-cm- 1  component. 

DISCUSSION 

Three  distinct  types  of  behavior  have  now  been  delin¬ 
eated  in  regard  to  the  intensity  variations  of  the  490-  and 
604-cm  ~ 1  peaks  from  fused  silica  as  follows:  (1)  Intensity 
increases  at  both  the  490-  and  604-cm  ~ '  positions  with  in¬ 
creasing  Tf,a  (2)  a  large  intensity  increase  at  604  cm  ~ 1  rela¬ 
tive  to  the  490-cm” 1  peak  observed  when  fused  silica  is  den- 
sified  by  neutron  irradiation,2  and  (3)  the  reversible  Raman 
intensity  increase  at  490  cm  ~ 1  observed  here.  It  should  be 
emphasized,  however,  that  the  type  1  and  2  variations  differ 
markedly  from  the  type  3  variation — annealing  is  required 
to  produce  reversibility  for  types  1  and  2. 

In  previous  Raman  work2"*  lines  near  490  and  604 
cm" 1  were  assigned  to  defects  of  the  type  — Si — O—Si — , 
where  the  dots  refer  to  a  broken  bond.  However,  in  view  of 
the  electron  paramagnetic  resonance  inactivity  of  ordinary 
bulk  fused  silica,17  it  is  more  appropriate  to  think  of  the  de¬ 
fect  as  a  stretched  bond,  or  related  to  the  effects  of  a 
stretched  bond.  Such  a  stretched  bond  would  have  a  smaller 
force  constant,  i.e.,  would  be  dynamically  weakened  but 
would  certainly  not  be  broken  from  an  energetic  point  of 
view,  that  is,  it  would  not  be  a  point  defect  in  the  usual  sense. 
Such  a  “defect"  could  be  considered  as  a  Si — O  bond  associ¬ 
ated  with  a  small  (e.g.,  less  than  120*)  bridging  bond  angle, 
and  consequently  with  a  large  (e.g.,  greater  than  1 .68  A)  Si — 
O  bond  length.'* 

In  regard  to  a  decreased  force  constant,  no  significant 
frequency  shifts  were  observed  in  this  work  and  in  previous 
Raman  work4  for  the  defect  lines  at  490  and  604  cm'1.  Also 
the  shifts  reported  for  the  defect  lines  by  Bates  el  at. 2  are  2% 
or  less.  Thus,  the  490-  and  604-cm  ~ 1  lines  probably  refer  to 
vibrational  modes  affected  by,  but  not  directly  related  to,  the 
Si — O  elongation  involved.  Further,  the  490-  and  604-cm  “ 1 
Raman  lines  are  unusually  sharp  for  glasses.  Thus,  it  is  not 
unreasonable  to  relate  these  lines  to  local  modes  decoupled 
from  the  silica  network  by  virtue  of  the  stretched  Si — O 
bond,  e.g.,  modes  of  the  03Si  or  OSiOj  units  and/or  of  the 
03Si — 0Si03  grouping,  in  which  one  intervening  and  elon¬ 
gated  Si — O  bond  replaces  the  dots  in  the  previous  — Si — 
O—Si — picture.  Because  the  concentration  of  stretched  Si — 
O  bonds  may  be  only  about  6  X  10*  Vcm3,4- ,v  or  about  0.1%, 
a  significant  lowering  of  the  stretching  frequency  at  1060 
cm  _  1  would  not  be  expected,  that  is,  most  of  the  Si — O  vi¬ 
brations  correspond  to  the  main  glass  structure. 

In  previous  Raman  studies  of  fused  silica,  it  was  esti¬ 
mated  that  a  small  fraction,  roughly  6x  10'Vcm'  of  de¬ 
fects,4  for  example,  of  highly  elongated  Si — O  bonds  might 
exist.  This  high-energy  fraction  would  be  "frozen”  into  the 
thermodynamically  unstable  structure  that  is  produced  by 
rapid  quenching.  Because  such  elongated  bonds  would  con¬ 
stitute  high-energy  sites,  they  would  probably  react  prefer¬ 
entially  with  water  to  form  adjacent  Si — OH  groups.”  Fur¬ 
ther,  high  tensile  stress  might  also  increase  the  fraction  of 


elongated  Si — O  bonds.  Or,  high  tensile  stress  could  alter 
structures  immediately  associated  with  the  elongated  bonds, 
i.e.,  the  effects  of  temperature  rise  and  of  tensile  stress  could 
be  very  different  mechanistically.  At  any  rate,  the  defects 
envisioned  here  are  topological  network  defects  associated 
with  elongated  high-energy  Si — O  bonds  whose  concentra¬ 
tions  may  be  changed  by  changes  in  TV.  or  alternatively 
through  preferential  chemical  reaction  with  H20  (Refs.  3 
and  4),  Ge02,  or  B203,3  etc. 

It  should  also  be  mentioned  that  when  uniaxial  stress  is 
applied  to  fused  silica,  changes  such  as  the  intensification  of 
Raman  shoulders  at  350-375  cm  ~ 1  and  at  1 1 5  ±  5  cm  - 1 
occur,  along  with  the  intensity  increases  at  490  cm' '.  Ap¬ 
parently,  the  main  silica  network  suffers  reversible  distor¬ 
tions,  as  well  as  the  reversible  bond  elongations  related  to  the 
defects.  The  intensifications  of  Raman  shoulders  observed 
here  are  significant,  and  do  not  appear  to  have  been  observed 
previously. 

The  magnitude  of  the  intensity  changes  observed  in  the 
tensile  stress  experiments  can  be  shown  to  be  surprisingly 
large  compared  to  the  energy  imparted  to  the  fused  silica 
optical  fiber  by  virtue  of  the  work  done  in  stretching  it.  Our 
measurements  indicate  that  a  7-kg  weight,  when  suspended 
from  a  fiber  200 /zm  in  diameter  (tensile  stress,  2. 19  GPa)  and 
447  cm  long,  stretches  the  fiber  about  12.7  cm.  The  work 
done  is  thus  0.41  keal/mol  Si02,  or  14.8  cal/cm3  Si02  This 
work  amounts  to  only  0.2%  of  the  standard  heat  of  forma¬ 
tion  of  Si02  glass  at  25  ”C,  (AH'}  =  202.5  keal/mol  SiOj),20 
and  thus  it  is  10"  3  times  smaller  than  the  Si — O  bond  energy. 
(The  Si — O  bond  energy  may  be  crudely  approximated  by 
dividing  the  AH}  by  4  because  there  are  4  Si — O  bonds  per 
stoichiometric  Si02.)  However,  if  the  defects  involve 
stretched,  rather  than  broken  bonds,  and  the  energy  due  to 
stretching  goes  disproportionately  into  the  defect  sites,  the 
effect  would  be  non-negligible.  In  this  regard,  the  following 
calculation,  although  forced,  is  instructive. 

Consider  that  the  defect  concentration  is  roughly 
6x  lO'Vcm3,4  and  that  this  concentration  increases  60- 
70%  due  to  tensile  stress,  as  might  be  inferred  from  the  in¬ 
crease  in  the  /<<w//440  ratio  (Table  II).  Then  the  change  in  the 
defect  concentration  would  be  about  4  x  10'Vcm3.  Because 
there  are  4  Si — O  bonds  per  stoichiometric  Si02,  the  total 
concentration  of  Si — O  bonds  is  about  9  X  10*l/cm3,  which 
means  that  the  fraction  changed  by  stretching  would  be 
4  x  10  ~ 4.  The  product  of  the  work,  about  1 5  cal/cm3,  times 
4  x  10~4  is  0.01  cal/cm3,  which  would  be  the  energy  impart¬ 
ed  by  equipartition  to  the  new  defects  contained  in  1  cm3. 
Because  this  amount  of  energy  seems  negligibly  small,  and 
the  Raman  intensity  change  at  490 cm" 1  is  visually  obvious, 
a  disporportionate  amount  of  the  work  of  stretching  might 
go  into  creating  new  defects.  The  remaining  energy,  which 
certainly  is  most  of  the  work  of  stretching,  would  go  into  the 
main  silica  network,  as  evidence  by  shape  changes  near  350- 
375  and  115  cm- '.  Alternatively,  because  TF  is  constant 
during  stretching,  a  small  fraction  of  the  work  of  stretching 
could  go  into  stretching  the  existing  temperature-induced 
defects  even  farther,  as  opposed  to  creating  new  defects.  This 
additional  stretching  could  produce  structural  changes 
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which  result  in  matrix  element  effects  leading  to  increased 
intensity  at  490  cm~ Unfortunately,  this  alternative  pro¬ 
cess,  although  almost  certainly  more  realistic  than  one  in¬ 
volving  production  of  additional  new  defects,  is  virtually  im¬ 
possible  to  approximate  energetically  at  this  time. 

With  regard  to  mechanistic  differences  between  Active 
temperature  and  uniaxial  tensile  stress,  the  reversible  nature 
of  the  present  experiments  should  be  strongly  emphasized. 
For  example,  annealing  is  required  to  remove  the  neutron- 
induced  defects  which  give  rise  to  greatly  increased  Raman 
intensity  at  604  cm-1.2  Whereas,  with  application  of  tensile 
stress,  the  Raman  spectral  changes  are  reversible  within  20- 
30  min,16  (although  the  actual  relaxation  times,  which  are 
probably  very  short  compared  to  20  min,  have  not  been  mea¬ 
sured16).  Further,  the  present  relaxational  behavior  may  also 
be  contrasted  to  the  situation  in  which  qualitatively  similar 
Raman  intensity  changes  in  the  490-  and  604-cm- 1  defect 
components  resulted  for  samples  having  TF  values  of 
1400  *C,  after  annealing  for  200  h  at  1 100  *C.4  Hence,  the 
presently  observed  reversibility  is  in  accord  with  a  bond 
stretching  mechanism,  whereas  bond  breaking  or  the  pro¬ 
duction  of  new  defects  would  probably  be  irreversible  at 
room  temperature. 

Since  completion  of  this  work,  Raman  experiments 
have  been  conducted  of  the  rate  of  increase  of  the  OH- 
stretching  peak  intensity  at  3700  cm'1,  relative  to  the  800- 
cm“ '  silica  peak  intensity,  for  a  silicone  rubber  clad  fused 
silica  optical  fiber,  with  and  without  tensile  stress.  To  stress¬ 
es  of  1 .56  GPa,  no  significant  increase  in  the  rate  of  the  OH- 
stretching  intensity,  above  the  slow  rise  characteristic  of  the 
water  uptake  of  the  unstressed  fiber,  could  be  detected. 
Thus,  the  reversible  increase  in  the  490-cm  ~ 1  peak  intensity 
observed  in  this  work,  and  the  lack  of  a  corresponding  effect 
on  the  rate  of  water  uptake,  contrasts  with  the  decreased 
intensity  observed  at  490  and  604  cm  " 1  when  the  OH  con¬ 
tent  increases  at  constant  Tf  *  It  would  appear  therefore  that 
no  additional  chemically  active  sites  result  when  increased 
tensile  stress  is  applied  to  a  fused  silica  fiber.  On  the  con¬ 
trary,  it  now  seems  that  the  490-cm  ~ 1  intensification  is  a 
matrix  element  effect  associated  with  a  structural  change  in 
existing  defects.  The  Raman  experiments  involving  the  rate 
of  water  uptake  will  be  described  more  fully  elsewhere. 

SUMMARY  AND  CONCLUSIONS 

The  results  of  a  new  type  of  experiment  are  presented  in 
which  high  tensile  stress  is  applied  to  a  long  fused  silica  opti¬ 
cal  fiber  while  simultaneously  obtaining  its  Raman  spectrum 
by  forward  scattering.  Intensification  of  a  sharp  Raman 
peak  near  490  cm  1  was  observed  with  increasing  tensile 
stress,  and  the  reversible  nature  of  the  intensification  con¬ 
trasts  with  previous  Raman  observations  in  which  annealing 
was  required  to  produce  reversibility.  The  Raman  intensifi¬ 
cation  was  interpreted  in  terms  of  a  decoupled  local  mode 
resulting  from  a  stretched  Si — O  defect  site. 
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ABSTRACT 

This  paper  illustrates  the  relationship  of  microstructure, 
composition,  and  ferroelectric  state,  to  the  crack  growth  resis¬ 
tance  of  ferroelectric-ferroelastic  materials  such  as  BaTi03  and 
Pb(Zr,Ti)03  (PZT).  Results  of  fracture  energy  determinations  are 
analyzed  based  on  predictions  and  observations  of  microcracking 
and  twinning.  Changes  in  the  ratio  of  Zr  to  Ti  in  PZT  for  compo¬ 
sitions  near  the  morphotropic  boundary  between  the  tetragonal  and 
rhombohedral  symmetries  are  shown  to  have  a  significant  effect  on 
fracture  energy.  Fracture  energy  data  in  PZT  are  interpreted  in 
terms  of  energy  absorbed  by  a  stress-induced  phase  transformation 
and  by  microcracking. 

INTRODUCTION 

At  grain  sizes  greater  than  «5  pm,  the  fracture  toughness 
of  both  BaTi03  and  PZT  has  been  shown  to  be  significantly  larger 
in  the  ferroelectric  state  than  it  is  in  the  paraelectric  state 
[1,2].  In  BaTi03,  this  toughening  has  been  attributed  to  the 
interaction  of  the  crack  with  ferroelastic  domains,  i.e. , 
twins. 
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Variations  of  the  fracture  toughness  of  PZT  compositions  near 
the  morphotropic  boundary  between  the  rhombohedral  and  tetragonal 
phases  has  been  attributed  to  effects  of  both  phase  transformations 
at  the  crack  tip  and  microcracking  [2].  Toughening  of  a  95*5 
PbZr03-PbTi03  composition  has  also  been  hypothesized  to  be  due  to  a 
combination  of  a  pressure- Induced  phase  transformation  and  micro¬ 
cracking  [3].  Thus,  the  observed  toughening  in  PZT  and  BaT10s 
comes  about  through  a  number  of  mechanisms.  The  purpose  of  this 
paper  is  to  summarize  the  available  information  on  the  effects  of 
microstructure  and  composition  on  the  strength  and  fracture  tough¬ 
ness  of  BaTi03  and  PZT,  and  to  try  and  assess  the  contribution  of 
each  possible  toughening  mechanism. 

EXPERIMENTAL  PROCEDURE 

The  materials  employed  in  this  study  were  either  obtained  com¬ 
mercially  or  prepared  by  relatively  standard  techniques  [1,2].  The 
critical  fracture  toughness  measurements  were  made  on  center- 
grooved,  double  cantilever  beam  specimens  «38  mm  x  12.5  mm  x  2  mm 
in  size.  An  ~  0.5  mm  saw  notch  served  as  the  starting  crack.  The 
specimens  were  loaded  to  failure  by  an  applied  moment  [4]  while 
immersed  in  silicone  oil  at  the  test  temperature  (25  °C,  ferroelec¬ 
tric  state,  150  °C,  paraelectric  state).  The  crosshead  rate  on  the 
test  machine  was  0.5  cm/min.  The  strength  measurements  were  made 
in  3-point  flexure  at  a  loading  rate  of  0.05  cm/min. 

RESULTS  AND  DISCUSSION 


In  BaTi03  the  critical  fracture  energy,  y,  is  a  function  of 
grain  size,  G,  in  the  ferroelectric  state  but  is  independent  of 
grain  size  in  the  paraelectric  state  (T^,  =  130  °C)  (Fig.  1).  The 

150  °C  data  fit  an  empirical  expression: 

y  =  -  0.02  Log  G  +  3.46  (1) 

The  small  value  of  the  slope  indicates  that  y  is  not  a  function  of 
grain  size  but  remains  constant  at  «3.5  J/m2.  Because  the  crystal 
structure  in  the  paraelectric  state  is  cubic,  this  type  of  behavior 
is  to  be  expected  [5].  The  only  deviation  to  this  trend  is  the 
data  obtained  on  a  Matopaammmrt*  two  phase  material  which  consists 
of  BaTi03  and  Mn02  (grain  size  =  30  to  50  pm).  A  y  of  12.8  J/m2 
was  measured  at  room  temperature  compared  to  a  value  of  12.7  J/m2 
measured  at  150  °C  on  this  material.  Fracture  analysis  of  this  two 
phase  material  showed  the  largest  amount  of  intergranular  failure 
of  any  of  the  BaTi03  materials  of  comparable  grain  size.  This 
difference  in  factor  morphology  suggests  a  different  controlling 
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size.  —  best  fit  to  room  temperature  data; 


...  contribution  due  to  twinning;  -  combined 

effect  of  twinning  and  microcracking  predicted  by 
Eq.  2. 

mechanism  for  crack  growth  in  this  two  phase  material  than  In 
single  phase  BaTi03. 

The  room  temperature  fracture  energy  of  BaTiOa,  having  a 
tetragonal  crystal  structure,  remains  constant  up  to  a  grain  size 
of  15  pm,  increases  rather  sharply  then  decreases  (Fig.  1).  This 
behavior  can  be  explained  on  the  basis  of  a  combination  of  twinning 
and  microcracking  enhanced  toughness  over  different  portions  of  the 
grain  size  range  as  follows. 

1  to  5  pm 

In  the  1  to  5  pm  grain  size  regime,  ferroelastic  twinning 
does  not  occur  [6],  so  the  fracture  energy  of  cubic  and  tetragonal 
BaTiOa  would  be  expected  to  be  quite  similar,  as  is  observed. 

These  are  the  lowest  values  of  y»  measured  experimentally. 

5  to  15  pm 

Above  a  grain  size  of  5  pm,  90°,  ferroelastic,  twins  are  pro- 
duced  by  the  cubic-to-tetragonal  phase  transformation  [7].  Based 
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on  the  increase  in  the  fracture  energy  of  single  crystal  BaTi03, 
from  0.8  to  1.4  J/m2  due  to  the  interaction  of  the  crack  with  the 
twins,  one  would  expect  the  y  of  the  polycrystalline  materials  to 
be  increased  proportionally.  Increases  of  this  magnitude  are  not 
observed  in  the  measured  values  of  y.  However,  it  is  observed  that 
fracture  over  this  grain  size  range  is  predominantly  intergranular. 
This  means  that  the  cracks,  instead  of  passing  through  the  twin 
structure,  take  the  easier  path  around  the  grains.  Since  no  crack- 
twin  interactions  take  place,  twinning  does  not  influence  the 
toughness. 

It  can  be  shown,  however,  that  microcracking  near  the  crack 
tip  should  affect  the  measured  fracture  energy.  That  microcracking 
must  play  some  role  in  the  fracture  process  is  shown  by  the  spon¬ 
taneous  cracking  of  polycrystalline  BaTi03  at  grain  sizes  between 
150  to  250  pm  [7].  Following  Rice  and  Freiman  [8],  one  can  calcu¬ 
late  the  contribution  that  microcracking  would  make  to  the  measured 
fracture  energy.  The  measured  y  can  be  calculated  from 

y  =  Tp  (1  -  §  )  +  H  [9EyBG)V2  -  Ae  EG]  (2) 


where  yp  is  the  fracture  energy  due  to  all  contributions  other 

than  microcracking,  i.e. ,  intrinsic,  twinning,  etc.;  G  is  the 
grain  size;  G  is  the  grain  size  at  which  spontaneous  failure 

occurs;  Ac  is  the  average  strain  mismatch;  yg  is  the  grain  bound¬ 
ary  fracture  energy;  E  is  Young's  modulus;  M  is  an  empirical  con¬ 
stant  related  to  the  number  of  microcracks.  Take  values  of  these 
parameters  for  BaTi03,  i.e.  E  =  12.5  x  lO10  Pa;  yg  =  single 

crystal  fracture  energy,  0.8  Pa-m;  At  =  5.6  x  10  4  based  on  the 

reduction  of  the  higher  local  internal  strains  due  to  twinning  [7]; 

1  4 

yp  =  3.5  Pa-m  for  G  <  «  18  pm,  and  =6.2  Pa-m  (3.5  x  gyg)  for 

G  >  w  18  pm;  Gg  =  200  pm,  and  M  =  2.5  based  on  the  best  fit  to 

fracture  energy  data  for  a  number  of  materials  [9].  The  y  calcu¬ 
lated  from  Equation  1  is  shown  as  the  dashed  line  in  Fig.  1.  Con¬ 
sidering  the  number  of  assumptions,  and  the  experimental  scatter  in 
both  the  measured  y  and  the  parameters  needed  to  calculate  it,  the 
fit  is  quite  good. 

15  to  50  pm 

Over  this  range,  both  microcracking  and  twinning  contribute  to 
the  measured  fracture  energy.  The  maximum  in  the  y-grain  size 
curve  can  be  calculated  from  the  following  expression  which  is 


obtained  by  differentiating  Equation  2  with  respect  to  G,  and 
setting  |^j  =  0. 


G 


max 


4  (y^ _  +  ZeE)2 

MZe  Gs 


(3) 


Using  the  same  values  of  the  parameters  as  before,  a  G  = 

MAX 

40  pm  is  calculated,  in  good  agreement  with  the  experimental  value. 
Note  that  this  calculation  Is  not  dependent  on  an  empirical  con- 
stant.  The  interaction  of  the  twins  with  the  propagating  crack  in 
a  material  having  a  grain  size  near  this  maximum  is  illustrated  in 
Fig.  2. 


50  to  200  pm 

At  grain  sizes  greater  than  about  50  pm,  both  twinning  and 
microcracking  still  contribute  to  increased  toughness,  but  linking 
of  microcracks  begins  to  occur,  leading  to  decreasing  fracture 
energy  with  grain  size.  The  fact  that  the  strength  of  BaT103 
having  an  average  grain  size  of  » 100  pm  is  not  lower  than  that  at 
fine  grain  sizes  [10]  does  not  rule  against  microcracking.  At  this 
grain  size,  flaws  are  contained  within  a  single  grain  and  propa¬ 
gate  at  the  single  crystal  fracture  energy.  Eventually,  spontan¬ 
eous  failure  will  take  place  leading  to  a  measured  y  that  is 
effectively  zero. 


Space  Charge  Effects 

Since  BaTi03  is  an  insulator  at  room  temperature,  it  was 
thought  that  space  charges  may  exist  in  the  material,  which  could 
create  additional  stresses  between  grains.  Semi-conducting  BaT103 
was  prepared  through  the  incorporation  of  samarium  in  solid  solu¬ 
tion.  These  materials  have  nearly  the  same  strains  due  to  the 
cubic-to-tetragonal  phase  transformations  but  cannot  maintain 
local  space  charges. 

The  fracture  energy  of  a  5  to  8  pm  grain  size  BaTi03  +  Sm  was 
4.1  J/m2  at  25  °C  and  4.5  J/m2  at  150  °C,  in  agreement  with  that 
determined  on  insulating  BaTi03.  Internal  stresses  calculated  from 
the  fracture  energy  and  the  flexural  strength  agreed  with  those 
determined  in  insulating  BaTi03  of  comparable  grain  size.  It  can 
be  concluded  that  any  space  charges  present  in  BaTi03  have  little 
or  no  effect  on  fracture. 
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PZT 


Table  1  summarizes  the  fracture  energy  and  strength  of  the 
PZT  compositions  investigated.  Note  that,  like  BaTi03,  the  frac- 
ture  energy  of  a  commercial  PZT  (TLZ-H)  is  significantly  higher  in 
the  ferroelectric  state  (25  °C)  than  in  the  paraelectric  (cubic) 
state  (200  °C).  As  with  BaTi03,  the  strength  of  this  material  is 
higher  in  the  paraelectric  state  due  to  the  presence  of  local 
tensile  stresses  at  25  °C  due  to  the  paraelectric-to-ferroelectric 
phase  transformation.  As  seen  in  Fig.  3,  the  strength  of  the 
material  in  the  two  states  decreases  and  crosses  with  increasing 
flaw  size  until,  at  some  flaw  size  large  enough  to  average  over 
the  compressive  and  tensile  components  of  the  internal  stress,  the 
ratio  of  the  two  strengths  is  essentially  equal  to  the  ratio  of 
1/2 

(y^/Vp)  as  would  be  predicted  by  fracture  mechanics  theory. 

For  purer,  experimental  materials,  there  is  a  distinct  minimum 
in  y  at  the  morphotropic  boundary  (53.5  PbZr03-46.5  PbTi03)  between 
the  tetragonal  and  rhombohedral  phases  (Fig.  4).  Less  extensive 
data  obtained  on  a  sintered  PZT  over  the  same  composition  range 
exhibit  similar  variations  in  y  (2). 

Table  1.  PZT  Fracture  Data 


Composition 

Pb(Zr.Ti)0, 

Crystal 

Structure 

(J/m2) 

Flexural  Strength 
MPa 

56 

Rhombohedral 

14.3 

163 

53.5/46.5 

Mixed 

6.2 

68 

50/50 

Tetragonal 

14.2 

104 

95/5 

Rhombohedral 

As  produced 

11.1 

58  [3] 

As  poled 

11.1 

60  [3] 

Pressure  Depoled 

17.7 

104  [3] 

** 

TLZ-H 

25  °C 

4.3 

81 

200  °C 

0.6 

96 

(paraelectric) 

** 

Plessey,  Australia 
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flaws  were  introduced  by  a  Knoop  indentor.  Ratio 
of  the  strengths  at  a  400  pm  flaw  size  is  propor- 

j- 

tional  to  the  ratio  of  (y)  for  the  two  conditions. 


The  variation  in  fracture  energy  with  composition  in  PZT  is 
interpreted  as  follows. 

At  the  morphotropic  boundary,  the  number  of  possible  twin 
states  is  maximized  because  states  for  both  the  tetragonal  and 
rhombohedral  symmetries  are  permitted.  Fracture  is  at  least  par¬ 
tially  transgranular.  The  interaction  of  the  crack  front  with  the 
twin  structure  within  the  grains  produces  increased  resistance  to 
propagation  in  analogy  with  the  processes  in  BaTi03  discussed  pre¬ 
viously.  As  seen  in  Fig.  5,  the  distortions  in  the  tetragonal  lat1 
tice  increase  with  distance  away  from  the  morphotropic  boundary. 
These  strains  would  be  expected  to  lead  to  microcracking.  Two 
pieces  of  evidence  for  this  phenomenon  are  the  extensive  cracking 
that  occurs  in  PbTi03  and  the  hysteresis  in  the  thermal  expansion 
curve  for  the  Pb(Tio.5,Zro.5)03  composition  as  reported  by  Biswas, 
et  al.  [11].  Similar  to  the  analysis  performed  on  BaTi03,  one  can 
calculate  the  contribution  that  such  microcracking  would  make  to 
the  fracture  energy  of  the  Pb(Ti0. 5 ,Zr0. s)03  composition.  Inser¬ 
ting  values  of  YpC  =6.2  Pa  m,  Ac  =  1.4  x  10"3  (measured  on  this 
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material),  G  =  200  pm,  G  =  8  pm,  E  =  81  x  10®Pa,  Yg  *  °-6  Pa-m,  and 

M  =  2.5,  into  Eq.  2,  one  obtains  a  value  of  y  58  9-4  J/"»2.  somewhat 
less  than  the  measured  value  of  14.2  J/m2.  For  a  constant  grain 
size,  one  can  also  derive  an  expression  which  predicts  the  composi¬ 
tion  in  the  tetragonal  phase  field  at  which  the  fracture  energy 
will  be  a  maximum.  If  one  assumes  that  the  elastic  modulus  is  con¬ 
stant  across  the  composition  range,  one  can  differentiate  Eq.  2 

with  respect  to  Ac.  By  setting  =  one  obtains  the  follow¬ 
ing  expression  '  ' 

in  ,)  *  ’•*  VEG>'/2  _  <4> 

For  the  values  of  Yg.  E»  and  G  used  previously,  4c^nax  ^  - 

1.4  x  10  3,  essentially  the  strain  expected  for  the  50/50  composi¬ 
tion.  y  1S  then  predicted  to  decrease  with  increasing  amounts  of 
T  i  O2  • 
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On  the  other  side  of  the  morphotropic  boundary  in  the  rhombo- 
hedral  phase  field,  the  lattice  distortions  are  much  smaller 
(Fig.  5).  This  is  borne  out  by  the  reduced  internal  stresses 
measured  in  the  56/44  composition.  Calculations  of  any  micro* 
cracking  contribution  to  y  for  this  composition  yields  a  smaller 
value  than  that  in  one  tetragonal  state  («8  J/m2).  Another  mech¬ 
anism,  namely  a  stress  induced  phase  transformation,  may  contribute 
to  toughening  in  this  composition.  It  is  expected  that  under  a 
tensile  stress  the  rhombohedral  phase  in  PZT  can  transform  to  the 
tetragonal  phase  [12].  When  this  transformation  occurs  at  the  crack 
tip,  it  can  reduce  the  strains  in  this  region  and  take  up  some  of 
the  energy  that  would  otherwise  have  gone  into  crack  growth.  How¬ 
ever,  no  quantitative  measure  of  the  contribution  of  this  mechanism 
of  toughening  can  be  given  at  present. 

Fracture  data  obtained  on  a  95  percent  PbZr03  -  5  percent 
PbTi03  composition  (  %  10  pm  grain  size)  also  suggest  a  contribu¬ 
tion  of  phase  transformations  to  toughening.  In  this  case  the 
material  was  tested  in  three  states,  as-produced,  poled,  and 
pressure  depoled.  The  material  that  was  depoled  by  the  application 
of  hydrostatic  compression  had  a  fracture  energy  that  was  about  50 
percent  greater  than  that  of  material  in  either  the  poled  or  as- 
produced  states  (Table  1).  The  fact  that  the  strength  of  the  95/5 


Figure  5.  Distortions  in  rhombohedral  and  tetragonal  lattices 
as  a  function  of  composition  in  PbZr03-PbTi03 
system  [7] 
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PZT  increases  after  pressure  deposing  suggests  that  phase  transfor¬ 
mations  play  a  more  important  role  than  microcracking  [3].  Further 
evidence  for  phase  transformation  induced  toughening  is  the  reduced 
fracture  energy  of  the  pressure  depoled  material  at  -60  °C.  No 
difference  in  y  was  measured  over  a  similar  temperature  range  for 
the  material  in  the  poled  or  virgin  states.  A  pressure- temperature 
phase  diagram  determined  by  Fritz  and  Keck  [13]  indicates  that  the 
transition  from  the  anti  ferroelectric  to  the  ferroelectric  state 
at  -60  °C  extends  into  the  tensile  stress  region.  This  shift  in 
the  ferroelectric-anti ferroelectric  boundary  suggests  that  more 
energy  would  be  needed  to  produce  a  phase  transformation  at  -60  °C 
than  at  room  temperature,  so  that  the  contribution  of  the  trans¬ 
formation  to  the  fracture  toughness  would  be  reduced. 

CONCLUSIONS 

Fracture  energy  data  for  both  BaTi03  and  PZT  have  been  inter¬ 
preted  in  terms  of  three  different  toughening  mechanisms.  Both 
twinning  and  microcracking  were  shown  to  contribute  to  the 
increased  room  temperature,  fracture  energy  of  BaTi03  for  grain 
sizes  greater  than  «15  pm.  The  grain  size  at  which  the  fracture 
energy  is  maximized  in  BaTi03  was  calculated  to  be  40  pm,  based  on 
a  model  of  Rice  and  Freiman,  in  good  agreement  with  observations. 

A  minimum  in  fracture  energy  was  measured  at  the  morphotropic 
boundary  between  the  rhombohedral  and  tetragonal  phases  of  PZT. 

It  was  hypothesized  that  while  twinning  was  the  major  contributor 
to  increased  fracture  energy  at  this  boundary,  the  significantly 
higher  y  in  the  tetragonal  phase  was  due  to  microcracking,  while 
the  equally  high  y  in  the  rhombohedral  phase  was  due  to  a  stress 
induced  phase  transformation.  Supporting  evidence  for  phase  trans¬ 
formation  toughening  in  a  95/5  PZT  was  also  given. 
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